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EXECUTIVE  SUMMARY 


Hie  success  of  achieving  a  Single  Integrated  Air  Picture  (SIAP)  is  a  function  of 
the  accuracy  with  which  the  positional  estimates  of  aerospace  objects  can  be 
made.  This  report  presents  a  detailed  analysis  of  how  track  accuracy  is  affected 
by  navigation,  sensor,  and  time  errors,  which  include  errors  due  to  time  syn¬ 
chronization,  time  stamping,  time  latency,  etc.  Although  the  analysis  assumes 
a  shipbome  sensor  system  for  illustrative  purposes,  the  analysis  is  general  and 
can  be  applied  to  all  sensing  systems  that  comprise  the  distributed  system  that 
develops  and  maintains  the  SIAP. 

Navigation  errors  are  the  result  of  errors  in  the  accelerometers  and  gyroscopes 
forming  the  inertial  navigation  system.  The  focus  here  is  not  on  the  errors  in 
the  individual  inertial  instruments  but  on  their  performance  as  they  affect  the 
uncertainty  in  the  position  of  the  inertial  navigation  system  and  the  orientation 
of  the  navigation  reference  frame.  Navigation  positional  uncertainty  or  error  is 
described  by  error  in  the  north  position,  east  position,  and  height  of  the  naviga¬ 
tion  system.  Error  in  the  orientation  of  the  navigation  reference  frame  is  de¬ 
scribed  by  three  misalignments  that  represent  the  cumulative  effect  of  the  gy¬ 
roscopes’  inability  to  maintain  a  perfect  inertial  reference.  Errors  in  navigation 
position  translate  directly  into  errors  in  track  measurement.  A  300  m  east  posi¬ 
tion  error  results  in  a  300  m  track  error.  Navigation  misalignments  also  have  a 
significant  effect  on  track  accuracy.  A  1  mrad  misalignment  about  all  three 
navigation  reference  axes  results  in  a  track  error  of  260  m  for  objects  100  nm 
away.  The  misalignment  contribution  is  proportional  to  range  to  the  object. 
Thus,  for  an  object  10  ran  from  the  ship,  the  1  mrad  misalignment  contribution 
is  26  m. 

Sensor  error  is  also  divided  into  two  groups  of  error.  First  are  errors  in  the  re¬ 
ported  location  of  the  object  relative  to  the  sensor  location.  This  error  is  quan¬ 
tified  by  error  in  the  reported  range  to  the  object  and  its  elevation  and  bearing 
relative  to  the  sensor’s  reference  axes.  The  second  group  of  sensor  errors  con¬ 
sists  of  those  relating  to  the  uncertainty  in  the  knowledge  of  the  sensor’s  refer¬ 
ence  axes.  Conceptually  these  errors  are  similar  to  navigation  misalignments  in 
that  they  represent  the  difference  between  the  computer’s  representation  of  the 
orientation  of  the  sensor  axes  and  their  true  orientation.  Sensor  range  error 
translates  directly  into  track  error.  Sensor  elevation  error  contribution  to  track 
error  is  proportional  to  the  range  to  the  object  while  the  bearing  contribution  is 
proportional  to  range  to  the  object  multiplied  by  the  cosine  of  elevation.  Track 
accuracy  sensitivity  to  sensor  misalignment  is  almost  identical  to  the  sensitiv¬ 
ity  to  navigation  misalignments  because  the  distance  from  the  sensor  to  the 
navigation  center  is  small  compared  to  the  range  to  the  object. 

Within  the  integrated  architecture,  potential  contributors  to  time  error  include 
clock  synchronization  error  (both  within  an  interfacing  unit  and  across  inter- 
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facing  units),  errors  in  the  time  stamping  of  measurements,  tracklets,  and 
track  messages,  processing  delays,  communication  network  (i.e..  Link  16)  de¬ 
lays,  and  uncertainty  in  the  knowledge  of  the  nominal  time  latency.  For  sim¬ 
plicity,  composite  time  errors  ranging  from  0.001  to  1  sec  are  considered.  The 
nominal  time  or  processing  latency  of  the  system,  which  is  defined  to  be  the 
nominal  time  the  track  information  is  used  after  it  is  observed  by  the  radar,  is 
shown  to  be  an  important  factor  for  objects  with  high  acceleration.  For  this 
analysis  a  range  of  potential  threats  is  considered  from  relatively  low-speed  ob¬ 
jects  (200  kt)  to  high-speed  combat  aircraft  (800  to  1,600  kt)  having  an  accel¬ 
eration  capability  on  the  order  of  9  g. 

The  sensitivity  of  track  error  to  time  errors  depends  on  the  object’s  velocity  and 
acceleration.  The  sensitivity  is  found  to  equal  the  vector  sum  of  the  velocity  and 
the  product  of  the  acceleration  and  the  time  latency  of  the  system.  For  a  high¬ 
speed  combat  aircraft  and  a  network  latency  of  5  sec,  the  time  error  is  magni¬ 
fied  by  a  factor  of  956  m/sec.  Thus,  a  time  error  of  1  sec  is  translated  into  a 
track  error  of  956  m.  The  time  latency  of  the  system  has  the  effect  of  increasing 
or  decreasing  the  contribution  of  the  object’s  acceleration  to  the  sensitivity.  For 
example,  reducing  the  time  latency  by  a  factor  of  two  can  reduce  the  track  er¬ 
ror  to  735  m. 

The  parametric  analysis  shows  that  worst  case  conditions  of  1  sec  total  com¬ 
posite  time  error  and  30  sec  network  latency  lead  to  track  errors  of  approxi¬ 
mately  3,000  m  for  high-speed  aircraft.  Composite  time  errors  on  the  order  of 
only  a  millisecond  lead  to  track  errors  of  less  than  4  m  for  these  worst  case 
conditions. 

A  nominal  track  error  budget  is  constructed  illustrating  the  error  contributions 
from  navigation,  sensor,  and  time.  This  budget  shows  a  balance  between  the 
navigation  and  sensor  contributions  to  track  error.  The  contribution  from  time 
error  is  a  distant  third  for  the  nominal  error  parameters  assumed.  The  time  er¬ 
ror  contribution  becomes  more  important  at  closer  ranges  as  the  other  range- 
dependent  contributions  decrease.  Sensitivity  curves  are  presented  which  show 
how  quickly  track  error  changes  as  a  function  of  each  of  the  individual  error 
terms. 

The  results  of  the  parametric  model  and  the  system  specific  error  parameters 
can  be  applied  to  using  and  testing  the  Integrated  Architecture  Behavior  Model 
(IABM)  within  the  Joint  Distributed  Engineering  Plant  (JDEP)  Technical 
Framework.  In  particular,  the  system  specific  error  parameters  define  the  sta¬ 
tistical  error  distribution  for  each  error  term.  The  parameter  values  coupled 
with  the  corresponding  distribution  define  the  level  of  bias  that  should  be  in¬ 
troduced  into  every  sensor  that  is  associated  with  an  instantiation  of  the  IABM. 
The  magnitude  of  the  contribution  of  each  error  to  track  accuracy  provides  a 
hierarchy  for  the  sequence  of  adding  errors  and  systems  into  the  IABM  for  test- 
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Should  be  introduced  into  the  disputed  system  beforo  “^Ztesys- 
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1. 


INTRODUCTION 


1.1  Problem  Statement 

A  critical  element  for  achieving  a  Single  Integrated  Air  Picture  (SIAP)  is  attain¬ 
ing  the  condition  of  correct  geodetic  alignment  of  positional  data  in  the  1984 
World  Geodetic  System  (WGS  84)  coordinate  frame.  Hie  degree  of  attaining  this 
geodetic  alignment  can  be  measured  by  the  accuracy  by  which  the  track  posi¬ 
tional  information  is  known. 

Track  accuracy  is  a  key  ingredient  in  how  well  units  in  a  distributed  system 
can  perform  their  various  functions.  For  example,  improvements  in  track  ac¬ 
curacy  improve  the  success  ratio  of  correlation  schemes,  which  in  turn  leads  to 
a  reduction  in  duals  and  a  reduction  in  misidentifications.  Reductions  in  track 
error  result  in  higher  values  of  track  quality  thus  improving  reporting  respon¬ 
sibility  throughout  the  distributed  system.  Improvements  in  track  accuracy 
and  its  subsequent  affect  on  improving  warfighter  effectiveness  in  a  distributed 
system  are  required  to  enable  advanced  concepts  of  netcentric  warfare  such  as 
cueing,  engage  on  remote,  and  post  before  processing. 

In  this  technical  report,  we  present  a  methodology  and  results  for  investigating 
and  understanding  track  accuracy.  In  particular,  the  methodology  provides  a 
straightforward  means  of  identifying  the  major  elements  that  contribute  to  er¬ 
rors  in  track  positional  estimates  and  in  quantifying  their  contribution  to  track 
error.  Because  precise  information  on  the  performance  of  sensors  is  unknown, 
we  present  the  results  of  a  parametric  study  that  illustrates  how  track  accu¬ 
racy  changes  as  the  errors  in  the  major  contributing  subsystems  change.  The 
particular  problem  we  address  is  that  of  a  ship  carrying  a  radar  sensor  that 
senses  an  aerospace  object,  Figure  1-1. 
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Aerospace  object 


Object  position  rela 
tive  to  earth  center 


Earth  center 


Object  position  rela¬ 
tive  to  sensor 


Object  position  relative  to 
ship’s  navigation  center 


Sensor  center 

Position  of  sen¬ 
sor  relative  to 
ship’s  navigation 
center 


Ship’s  navigation 
center 
Position  of  the  ship’s 
navigation  center  rela¬ 
tive  to  earth  center 


Figure  1-1  SIAP  Track  Accuracy  Problem 


1.2  Technical  Approach 

The  analysis  approach  is  to  develop  the  equation  that  maps  the  location  of  the 
aerospace  object  expressed  in  the  sensor  frame  to  the  location  of  the  object 
relative  to  the  Earth  center  expressed  in  the  WGS  84  coordinate  frame,  and  to 
consider  how  measurement  errors  and  misalignments  affect  the  estimate  of 
track  position. 

In  general,  sensor  measurements  are  made  relative  to  a  specified  coordinate 
frame.  For  example,  radar  sensor  measurements  are  made  relative  to  the  sen¬ 
sor  coordinate  frame,  the  ship’s  inertial  navigation  system  (INS)  provides 
measurements  of  position  and  orientation  relative  to  an  inertial  coordinate 
frame,  and  ship  structures  are  referenced  to  a  body  or  vehicle  frame  that  can 
rotate  relative  to  the  inertial  navigation  frame.  Meaningful  information  can  be 
exchanged  between  a  single  unit’s  various  sensors,  or  the  ship  and  other  units 
in  a  distributed  system  only  if  the  location  and  orientation  of  the  various  coor¬ 
dinate  frames  are  known  relative  to  each  other. 

An  examination  of  the  mapping  of  an  object’s  position,  expressed  in  the  sensor 
frame,  to  its  position  expressed  in  the  WGS  84  frame  identifies  three  groups  of 
error  contributors: 
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•  Navigation  position  (north,  east,  and  height)  and  misalignments 

•  Sensor  measurement  (range,  bearing,  and  elevation)  and  mis¬ 
alignments 

•  Time 

After  the  error  mechanisms  are  identified,  a  first-order  linear  perturbation 
analysis  provides  the  sensitivities  of  track  accuracy  to  the  identified  error  con¬ 
tributors.  Once  the  sensitivities  are  calculated,  the  error  budgets  and  paramet¬ 
ric  sensitivities  are  determined. 

The  approach  is  general  and  can  be  applied  equally  well  to  land-based  and  air¬ 
borne  sensor  systems.  The  analysis  presented  in  the  remainder  of  the  report  is 
intended  to  be  complete,  including  statements  of  assumptions  and  technical 
details.  An  examination  of  the  assumptions  provides  insight  into  how  the 
analysis  can  be  extended  to  other  factors  that  affect  track  accuracy. 

1.3  Motivation 

The  error  budget  approach  and  analysis  presented  in  this  report  is  a  powerful 
technique  for  characterizing  and  understanding  the  performance  of  a  family  of 
systems  such  as  the  distributed  architecture  representation  of  the  SIAP.  It  can 
be  applied  to  existing  systems  to  understand  what  drives  current  performance 
and  suggest  improvements,  and  it  can  be  applied  to  systems  in  the  design 
phase  to  predict  performance  and  perform  trade-off  studies. 

The  error  analysis  process  is  a  rigorous,  quantitative  examination  into  how  a 
family  of  systems  performs.  Its  power  comes  from  the  fact  that  it  can  be  used 
to  examine  the  total  family  of  systems,  as  well  as  individual  systems,  subsys¬ 
tems,  and  single  components.  The  process  involves  identifying  all  mechanisms 
that  introduce  error,  tracing  the  propagation  of  the  error  throughout  the  family 
of  systems,  and  quantifying  the  effect  on  system  performance. 

The  error  is  a  combination  of  two  parts:  the  deterministic  error,  which  can  be 
compensated,  and  the  random  error,  which  in  most  cases  can  be  reduced 
through  improvements  in  the  build  process  or  through  external  measurements 
and  corrections.  Error  budget  analysis  provides  an  assessment  of  the  perform¬ 
ance  of  the  family  of  systems  before  improvements  are  made.  It  can  identify 
candidate  improvement  options  and  provide  projections  of  performance  if  com¬ 
pensation  techniques  are  applied  or  improved  build  processes  are  initiated. 
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2.  POSITION  ERROR  ANALYSIS  FOR  AEROSPACE  OBJECTS 


2. 1  Introduction  to  Position  Error  Analysis 

Imagine  some  weapon  system  that  carries  a  sensor  system  that  senses  an 
aerospace  object— think  of  a  ship  that  carries  a  radar  that  senses  an  aircraft 
The  position  and  attitude  of  the  navigation  reference  frame  relative  to  the  Earth 
are  measured  by  an  Inertial  Navigation  System  (INS).  There  is  a  lever  arm  be¬ 
tween  the  center  of  the  INS  and  the  sensor  system,  for  they  may  be  located  a 
considerable  distance  apart  on  the  ship.  The  sensor  system  measures  the 
range  and  angles  to  the  object,  angles  with  respect  to  some  sensor  coordinate 
frame,  which,  in  turn,  are  related  to  some  ship  referenced  coordinate  frame. 
Given  the  sensor  system  errors,  the  sensor-to-ship  orientation  errors,  the  lever 
arm  errors,  and  the  INS  errors  in  both  position  and  attitude,  what  are  the  re¬ 
sulting  errors  in  absolute  object  location?  That  is,  what  are  the  errors  in  the 
location  of  the  object  with  respect  to  the  Earth,  not  just  with  respect  to  the 
ship? 

The  analysis  here  is  basic  and  generic.  It  does  not  apply,  in  detail,  to  any  spe¬ 
cific  set  of  operational  systems,  but  it  is  applicable,  in  general,  to  many  such 
combinations.  The  objective  here  is  a  set  of  definitions  and  a  consistent  nota¬ 
tion  that  together  lead  to  an  understanding  of  error  propagation  and  to  equa¬ 
tions  useful  for  error  estimation  and  for  simulation. 

In  the  operational  world,  the  process  begins  with  the  sensor  measurements  of 
an  object  to  relate  the  object  s  location  to  the  sensor,  then  the  object’s  location 
is  related  to  the  ship  s  INS,  and  finally,  the  object’s  location  is  related  to  the 
Earth.  The  discussion  of  the  process  begins  in  reverse,  starting  with  the  Earth, 
then  turning  to  the  ship’s  INS,  then  to  the  sensor,  and  finally  to  the  object.  In 
the  development  of  the  results,  the  definitions  and  the  error  relationships  are 
introduced  as  they  are  needed.  The  equations  are  in  a  logical  order  for  devel¬ 
opment  and  discussion,  but  not  necessarily  for  use.  Summary  sections  are  in¬ 
cluded  to  compile  the  definitions  and  the  measurement  processing  steps  into  a 
more  compact  reference  form  that  may  assist  in  understanding,  and  to  restate 
the  equations  in  a  more  usable  way.  Finally,  the  equations  are  reordered  again, 
with  additions,  to  outline  a  sequence  of  computations  for  use  in  simulation  and 
analysis. 


2-1 


7.2.6.3_TR(4-002)_1.0„JSSEO_040407 


2.2  The  Earth 


Consider  an  ellipsoid  model  of  the  Earth  geoid  L  in  particular,  the  one  used  as 
geodetic  datum  by  the  World  Geodetic  System  1984  (WGS  84).  Parameters  re¬ 
lated  to  this  ellipsoid  are  listed  in  Table  2-1. 


Table  2-1  WGS  84  Ellipsoid  Parameters 


Parameter 

Notation 

Formula 

Value 

Angular  Velocity 

<»e 

7.2921151467 

xl0~5rad/sec 

15.04106718  dee/hr 

Semimajor  Axis 

ae 

6378137  m 

Flattening 

f 

1/298.2572236 

Semiminor  Axis 

K 

6356752.3  m 

Eccentricity 

e 

V/(2-/) 

0.08181919084 

Eccentricity 

Squared 

e 2 

e2 

0.00669437990 

Axis  Ratio 

bjae 

yll-e2 

0.9966471893 

Axis  Ratio  Squared 

1-e2 

0.9933056199 

Linear  Eccentricity 

E 

eUe 

521854.01  m 

Minor  Eccentricity 

e 

0.08209443796 

Let  the  geodetic  coordinates  of  a  point  B  with  respect  to  this  datum  be  denoted 
as  (<p,A,h) ,  which  are  the  geodetic  latitude,  the  longitude,  and  the  altitude  of  B 

above  the  ellipsoid.  The  geodetic  latitude  and  altitude  are  illustrated  in  Figure 
2-1.  6 


The  WGS  84  geoid  is  that  particular  equipotential  surface  of  the  earth  that  coincides  with 
mean  sea  level  over  the  oceans  and  extends  hypothetically  beneath  all  land  surfaces.  The  geoid 
serves  in  land  areas  (if  conventional  leveling  is  not  available)  as  the  vertical  datum  or  reference 
surface  to  which  height-above-mean  sea  level  values  are  referenced.  In  ocean  areas  the  geoid 
serves  as  the  geodetic  height  of  points  on  the  ocean's  surface. 
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polar  axis  ellipsoid 


Figure  2- 1  Geodetic  Latitude  and  Altitude 

Define  the  usual  Cartesian  coordinate  system,  called  the  e-frame,  or  the  Earth 
frame,  to  be  the  Conventional  Terrestrial  System  (CTS)  associated  with  the 
WGS  84  datum.  (All  Cartesian  coordinate  systems  used  here  are  orthogonal, 
right-handed  systems.)  This  frame  is  fixed  to  the  ellipsoid,  and  so  to  the  Earth. 
Its  center  is  at  the  ellipsoid  center  (called  here  the  center  of  the  Earth,  or  the 
point  O).  Its  z-axas  is  along  the  ellipsoid  axis  of  symmetry  (called  here  the 
Earth  rotation  axis,  or  north  pole,  which  ignores  polar  motion).  Its  x-axis  and 
y-axis  are  in  the  plane  of  ellipsoid  symmetry  known  as  the  equatorial  plane. 
The  x-axis  is  in  the  Greenwich  meridian  (the  meridian  of  zero  longitude).  Earth 
coordinates  are  depicted  in  Figure  2-2. 

Coordinate  systems  here  are  assigned  definite  centers  for  ease  of  conceptuali¬ 
zation,  but  they  are  not  used  as  systems  of  positional  coordinates.  Instead, 
they  are  used  as  systems  of  vector  coordinates,  so  their  centers  do  not  enter 
into  transformations  from  one  system  to  another.  It  is  only  the  directions  of  the 
coordinate  axes  that  matter.  That  is,  all  the  vectors  of  “position”  here  are,  in 
fact,  displacement  vectors  that  go  from  one  specified  point  in  space  to  another. 
Vectors  conceptualized  without  regard  to  coordinate  systems  are  called  geomet¬ 
ric  vectors. 
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A 


Figure  2-2  Earth  Coordinates:  the  e-Frame 


A  geometric  vector  is  never  equated  to  its  coordinate  system  representation  as 
a  mathematical  3-vector.  That  relationship  is  expressed  by  the  notation 
r  *-»  (r)e ,  which  means  that  the  geometric  vector  r  is  expressed  by  or  is  repre¬ 
sented  by  the  coordinate  vector  (r)e  in  the  e-frame.  This  approach  adds  some 

notational  complexity,  but  it  ultimately  serves  to  reduce  confusion  when  sev¬ 
eral  coordinate  systems  are  in  use  simultaneously,  as  they  are  here. 

Let  rh  be  the  vector  from  the  center  of  Earth  to  B  .  The  components  of  rb  in  the 
e-frame  are  given  by 


( = 


(r™  +  h)CfCx 
(rm,  +  h)CfSx 


(2.2-1) 


where  C,  and S#  are  shorthand  for  co$<p  and  sin  <f> ,  etc.;  and  rm,  is  the  east-west 
radius  of  ellipsoid  curvature,  a  function  of  the  latitude,  given  by 

\-eXrm  (2.2-2) 


2.3  The  Ship 

Consider  a  ship,  fixed  or  moving,  in  the  vicinity  of  the  Earth.  Define  a  point  B  , 

fixed  to  that  ship,  called  the  center  of  the  ship’s  INS,  with  geodetic  coordinates 
as  above. 
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The  geodetic  coordinates  of  the  INS  are  provided  by  the  inertial  instruments 
forming  the  INS,  and  are  subject  to  error.  Let  W,X,h)  be  the  measured  geodetic 
coordinates  of  B ,  and  define  the  geodetic  coordinate  errors  (S0,SX,8h) ,  so  that 
the  measured  quantities  are  the  true  quantities  plus  the  errors,  or 

<j>  =  $  +  8<p ,  X-X  +  SX,  h=h  +  8h  (2.3-1) 

Measured  values  here  are  denoted  with  a  tilde  R;  true  values  are  “plain”  sym¬ 
bols,  and  error  values  have  a  preceding  8 . 

The  measured  coordinates  of  (rh  )e ,  in  terms  of  the  true  values  and  the  errors, 
are 


(Fbh=(rb)e+(S?b)e  (2.3-2) 

where,  approximately  (to  first  order  in  the  8  quantities) 


+  h)8<p 

-V, 

+  (rew  +  h)C^SX 

f 

L  . 

+  8h 

c#c/ 

e 

- 0  J 

e 

S.  J 

where  is  the  north-south  radius  of  ellipsoid  curvature,  a  function  of  the  lati¬ 
tude,  given  by 

rns  =  a({\~e2){\-e2SlYV2  =  rew{\-e2){\- e2Sl)~'  (2.3-4) 

Define  a  second  Cartesian  coordinate  system,  called  the  1-ffame,  or  the  local 
level  frame,  to  have  its  center  at  B ,  and  its  x-,  y-,  and  z-axes  along  the  north 
east,  and  down  directions  at  B  .  That  is,  the  axes  of  the  local  level  frame  have 
directions  along  the  unit  vectors  (N,  E,  D ) .  Down,  here,  means  orthogonal  to  the 

ellipsoid  at  the  surface  of  the  ellipsoid.  In  the  Earth  frame,  the  components  of 
the  unit  vectors  along  the  local  level  frame  axes  are 


A 

-V/ 

'-si 

-C& 

(N)e~ 

.  (E)e  = 

.  0)e  = 

-c.s. 

.  c. 

e 

0 

e 

9  Je 
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The  local  level  system  is  not  fixed  to  the  Earth  (unless  the  ship  is  fixed  to  the 
Earth),  and  it  is  not  fixed  to  the  ship;  its  axes  rotate  as  the  ship’s  position 
changes.  Figure  2-3  illustrates  the  local  level  frame. 


Figure  2-3  Local  Level  Coordinates:  the  1-Frame 

The  errors  in  the  measured  components  of  the  INS  position  vector  are 

(<%),  =  (r„  +  h)8ftN)t  +  (rw  +  h)C,5Z(E)e  - Sh(D)e  (2.3-6) 

That  is,  an  error  in  latitude  produces  an  error  in  the  northward  direction,  an 
error  in  longitude  produces  an  error  in  the  eastward  direction,  and  an  error  in 
altitude  produces  an  error  in  the  (negative)  downward  direction. 

The  rotation  matrix  that  transforms  vectors  from  expression  in  the  e-frame  to 
expression  in  the  1-frame  is  given  by 

(2.3-7) 

where  the  standard  trio  of  single-axis  rotation  matrices  is  defined  as 


Rk  =R2(- 


n 


-MM)-- 


s,cx 


cc\ 


-v, 

Cjl 

c,sx 


0 
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^  0  0  cos#  0  -sin#  cos#  sin#  0 

/?](#)=  0  cos#  sin#  ,  R2{0)=  o  1  0  ,  /?,(#)  =  -sin#  cos#  0  (2.3-8) 

0  -sin#  cos#  sin#  0  cos#  0  0  1 

Notice  that  the  rows  of  Rle  are  the  elements  of  the  three  local  level  frame  axis 
unit  vectors  expressed  in  the  Earth  frame. 

For  example,  the  INS  position  vector,  expressed  in  local  level  coordinates,  is 
given  by 

(rb),  =K(h)e  =  0  (2.3-9) 

This  vector  has  a  small  northward  component  (because  the  Earth  is  elliptical, 
not  spherical),  no  eastward  component,  and  a  large,  negative,  downward  com¬ 
ponent,  as  might  be  expected. 

The  errors  in  latitude  and  longitude  produce  errors  in  the  transformation  from 
the  Earth  frame  to  local  level  coordinates.  That  is 

Rle=Rle+SRle  (2.3-10) 

To  first  order  in  the  small  angle  quantities,  we  relate  the  error  in  the  transfor¬ 
mation  to  the  measured  transformation  by  the  antisymmetric  matrix,  5Qle , 
such  that 

5Rle=5QleRle  (2.3-11) 

and 

Snie  =-^Q2  +SJUt2(-?~-0)Q}Rt(-~-0)  (2.3-12) 

where  the  standard  trio  of  single-axis  antisymmetric  matrices  are  defined  as 

o  o  ol  To  o  -il  ro  i  o' 

Q,  =  0  0  1  .  n2  =  0  0  0  ,  Q3=  -1  0  0  (2.3-13) 

0  -1  0  1  0  0  0  0  0 

It  helps  in  the  derivation  of  the  result  above  to  note  that  to  first  order  in  56 

Rt  (#  +  56)  —  Rt  (6)Rj  (56)  =  R,  (#)  +  Ri  (0)C1.56 ,  Ri(6)Qi  =  Q.  /?.(#)  (2.3-14) 

In  total  then,  the  small  antisymmetric  matrix  5Qle  is  associated  with  the 
small  3-vector  (50)t ,  where 
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&p2 

=  -5$ 

'o' 

1 

+  SA 

© 

■ - 1 

_<^3_ 

i 

0_ 

i 

-s 

*  J/ 

(2.3-15) 


(It  makes  no  meaningful  difference  when  computing  small  error  quantities, 
such  as  the  one  above,  whether  the  true  or  the  measured  values  of  the  parame¬ 
ters  are  used.)  That  is 


SQle  (S$), , 


0 


&/>> 

0 


-&t>2 

0 


(2.3-16) 


(The  symbol  <=>  denotes  that  an  antisymmetric  matrix  and  a  3-vector  are  asso¬ 
ciated,  as  they  are  in  the  relationship  above.) 

A  warning  to  the  reader  is  in  order.  When  two  (or  more)  coordinate  systems  are 
involved,  say  the  a-frame  and  the  b-ffame,  and  there  is  some  orthonormal 
transformation  matrix,  say  Rha ,  that  transforms  vectors  from  expression  in  the 
a-frame  to  expression  in  the  b-frame,  then,  for  any  “whole”  vector,  say  v  ,  it  is 
true  that  its  two  coordinate  expression  3-vectors  satisfy  the  relationship 


(v)*=*L(v). 


(2.3-17) 


When  there  are  errors  in  the  coordinates  of  this  vector,  those  errors  are  typi¬ 
cally  denoted  as  (Sv)a ,  and  (Sv)h ,  but  it  is  not  generally  true  that 


(Sv)b  =  R," (< 5v)a  NOT  GENERALLY  TRUE  (2.3- 18) 

This  discrepancy  happens  because  there  are  two  sources  of  error  in  the  com¬ 
ponents.  One  source  is  the  error  in  the  initial  coordinate  system.  This  error 
does  transform  according  to  the  equation  above.  The  other  source  is  due  to  er¬ 
rors  in  the  transformation  itself,  which  adds  an  extra  term  to  (5v)b ,  so  that  the 
correct  relationship  is  usually 

(Sv)b  =  Rba(Sv)a+SRba(v)a  (2.3-19) 

That  is,  (Sv)a  and  (Sv)b  are  not  the  coordinates  of  the  same  vector,  although 

the  notation  makes  them  look  so.  The  notation  could  be  modified  to  accommo¬ 
date  this  issue,  but  the  result  is  so  ugly  and  confusing  itself  that  it  is  not  done 
here. 


Fixed  to  the  ship,  centered  at  the  ship’s  INS  center,  B ,  is  yet  another  Cartesian 
coordinate  system,  called  the  b-frame,  or  the  body  frame,  having  the  subscript 
b  to  denote  its  coordinate  vectors.  It  is  convenient  that  this  body  coordinate 
system  be  coincident  with  the  local  level  frame  under  some  “nominal”  condi- 
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tions.  For  a  moving  system  such  as  a  ship,  the  nominal  conditions  are  that  the 
ship  be  “level”  and  be  pointed  in  the  north  direction.  So  the  x-axis  would  point 
forward,  the  y-axis  would  be  starboard  (out  the  right  side),  and  the  z-axis 
would  point  toward  the  bottom  of  the  ship.  The  b-frame  is  illustrated  in  Figure 
2-4. 


Figure  2-4  Body  and  Sensor  Coordinate  Frames 

The  exact  relationship  of  the  body  frame  to  the  local  level  frame  is  provided  by 
three  Euler  angles.  (y/b,Ob,<j>h) ,  which  are  called  the  (true)  body  heading,  body 
pitch,  and  body  roll. 

The  rotation  matrix  that  transforms  vector  coordinates  from  the  1-frame  to  the 
b-frame  is 


II 

(2.3-20) 

'V 

I 

to 

1 

S*  +C„Cn 

(2.3-21) 

CACn  +  Vr, 

Like  the  geodetic  coordinates  of  the  body  frame,  the  body  Euler  angles  are 
measured  quantities,  presumably  provided  by  the  INS,  and  are  subject  to  error, 
so  the  measured  quantities,  the  true  quantities,  and  the  error  quantities  are 
related  by 
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(2.3-22) 


¥b~V/h+  $Yh  ’  0h~0h+S0h,  0k=<Pb+S0h,  Rhl  =Rhl  +  SRU 


It  follows  that,  to  first  order  in  the  small  angle  quantities,  there  is  some  anti¬ 
symmetric  matrix,  SO.,, ,  such  that 


SRh,  —  SQhlRhl 


(2.3-23) 


and 


SQbI  =  +  86hRt  (0h  )QX  (A )  +  Sy/hRx  (0b  )R2  (§b  )Q  3Rf  (A  )K  (A )  (2.3-24) 

In  total  then,  the  small  antisymmetric  matrix  8Qhl  is  associated  with  the  small 
3-vector  {80)  h  of  navigation  misalignments  or  tilt  errors,  where 


That  is 


(2.3-25) 


sr 

ti 

80 u 

S0E 

Y 

0 

+  80h 

P  O 

+8ysh 

'-5 

80 o _ 

b 

0 

b 

rs*l 

b 

(3QW  <=>  ( 86)h , 


0 

80, 


S0D 

0 


S0F  - 86 ; 


-Me 

S0N 

0 


(2.3-26) 
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2.4  The  Sensor  System 

Also  (nominally)  fixed  to  the  ship  is  some  sensor  system,  for  example,  a  radar. 
Define  the  point  S  to  be  the  “center”  of  the  sensor  system.  The  displacement 
vector  from  B  to  S  will  be  denoted  by  A7hs ,  and  is  called  the  body-sensor  lever 

arm  vector.  This  vector  is  assumed  measured  in  body  coordinates.  Its  errors 
will  be  considered  later. 

The  sensor  system  is  presumed  to  have  its  own  Cartesian  coordinate  system, 
called  the  s-ffame,  or  sensor  coordinates.  The  rotation  matrix  that  transforms 
vectors  from  expression  in  sensor  coordinates  to  expression  in  body  coordi¬ 
nates  is  denoted  as  Rbs .  Nominally,  this  rotation  matrix  is  constant,  because 

the  sensor  system  is  assumed  fixed  to  the  ship.  In  practice,  however,  some  sys¬ 
tems  are  subject  to  deformation.  Ships,  as  an  example,  deform  for  many  rea¬ 
sons,  including  the  rearrangement  of  fuel,  stores,  or  cargo;  temperature  differ¬ 
ences  caused  by  the  sun  or  engines;  and  stresses  caused  by  propulsion,  sea 
state,  wind,  current,  or  even  armament  usage.  So  the  relationship  between  the 
sensor  frame  and  the  system’s  body  frame  might  not  be  constant.  Constant  or 
not,  its  measurement  is  considered  to  be  prone  to  error. 

Figure  2-4,  above,  illustrates  the  s-ffame  and  the  body-sensor  lever  arm.  In  the 
diagram,  the  sensor  frame  and  the  body  frame  are  nominally  identical  (in  axes 
directions).  These  two  coordinate  frames  might  have  any  directional  relation¬ 
ship  to  one  another,  and  that  relationship  will  not,  itself,  be  parameterized  here 

by  Euler  angles,  just  by  the  direction  cosine  (transformation)  matrix  mentioned 
above. 

Denote  the  measured  transformation  matrix  from  sensor  to  body  coordinates, 

K  •  and  the  error  in  this  transformation,  5Rhi ,  so  that 


(2.4-1) 

It  must  be  that 

(2.4-2) 

where  we  have  defined  the  matrix 

<5£2»,  =a-*J?>=<wM£ 

(2.4-3) 

It  follows  that 

sa,a+sal^SRtMi 

(2.4-4) 
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so  SQhs  is  antisymmetric  to  first  order  in  the  error  quantities,  and  we  may  as¬ 
sociate  it  with  a  3 -vector,  as 


0  Sy/2  ~Sy/2~\  lSy/ 

-Syr 3  0  Syr ,  =  <=>  (£$!/%  =  (2.4-5) 

_  <V2  -Syr ,  0  J  fc 


It  is  the  3-vector  above  that  will  characterize  the  errors  in  the  relationship  of 
sensor  coordinates  to  body  coordinates.  These  errors  will  be  called  the  sensor 
misalignments. 
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2.5  The  Aerospace  Object 


Somewhere,  remote  from  the  ship  and  the  sensor  system  it  carries,  is  an  aero¬ 
space  object  that  is  sensed.  Let  T  be  a  point  fixed  to  the  object,  called  the  ob¬ 
ject  center,  presumably  the  very  point  sensed  by  the  sensor  system.  Let  the 
vector  Arsl  be  the  displacement  from  the  center  of  the  sensor  to  the  center  of 

the  object,  called  the  sensor-object  displacement  vector.  Assume  that  the  sen¬ 
sor  system  measures  the  range,  the  bearing,  and  the  elevation  of  the  object 
with  respect  to  sensor  coordinates.  Figure  2-5.  Name  (the  true  values  of)  these 
parameters  (p,y/s,6s). 


Figure  2-5  Sensor  Measurement  Items:  Range,  Bearing,  and  Elevation 

Then 

(AK,)S  =P 

Of  course,  these  sensor  parameters  are  measured  with  errors,  so,  as  usual  here 

p  =  p  +  Sp,  WS=W,  +Sy/s,  ds=0s+Ses  (2.5-2) 

which  produce  errors  in  the  sensor-object  displacement  vector  denoted  as 

(4),  +(#*).  (2.5-3) 

and  given  by 


C0,C¥, 


(2.5-1) 
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(3rJs=fy 


Ce,Cv' 

Xcr." 

'-c9s9; 

-pWs 

+  P$VS 

c,c¥t 

.  . 

s 

ce 

"  J 

s 

0 

(2.5-4) 


The  sensor  object  displacement  vector  is  transformed  from  sensor  coordinates 
into  body  coordinates  as 


=/?fa(ArJ,)i  (2.5-5) 

Because  there  are  errors  in  both  the  measured  sensor-object  displacement  vec¬ 
tor  in  sensor  coordinates  and  in  the  transformation  from  sensor  coordinates  to 
body  coordinates,  there  are  conceptually  two  (groups  of  three)  sources  of  errors 
in  the  measured  sensor-object  displacement  vector  in  body  coordinates.  That  is 

(A I, )„  =  (Ar„ ),  +  (Srst  )b  (2.5-6) 

where 

(dfs,  )b  =  Rbs  ( Srs,  )s  +  SRbx  (A rs,  )s  (2.5-7) 

Next,  in  the  chain  of  object  vector  processing,  the  lever  arm  from  the  body  (or 
INS)  center  to  the  sensor  system  center  is  added  to  the  sensor-object  displace¬ 
ment  vector  to  create  the  body-object  displacement  vector,  which  goes  from  the 
point  P  to  the  point  T  .  In  symbols 

(Ar„  )b  =  (A  rbs  )b  +  (A r„  )b  (2.5-8) 

That  is,  (A rhl  )b  is  the  newly  defined  body-object  displacement  vector,  expressed 
in  body  coordinates.  There  are  three  (group)  sources  of  error  in  the  body-object 
displacement  vector;  the  two  groups  that  contribute  to  the  errors  in  the  sensor- 
object  displacement  vector,  and  the  errors  in  the  knowledge  of  the  lever  arm. 

The  lever  arm  errors  are  described  directly  by 

(AFfo  )h  =  (A  rbs  )b  +  (Srbs  )b  (2.5-9) 

where,  of  course,  the  errors  themselves  are  represented  by  (SFhs)b.  In  symbols, 

(A?bI  )b  =  (A ?bt  )b  +  (Srbt  )h  (2.5-10) 

where  the  total  error  in  the  body-object  displacement  vector  is 

(K  =  (K  h  +  {Sr„  )b  (2.5-11) 

The  next  step  in  object  vector  processing  is  to  transform  the  body-object  dis¬ 
placement  vector  from  body  coordinates  to  local  level  coordinates,  as 

(A^,)/  =Rl,(Arbt)b  (2.5-12) 
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Errors  in  measured  body  attitude  (roll,  pitch,  and  yaw)  cause  errors  in  the 
transformation  matrix  from  local  level  to  body  coordinates,  and  thus  add  errors 
to  the  body-object  displacement  vector  in  local  level  coordinates.  In  symbols 

(AF*),=(A  rbt),+(Srbl),  (2.5-13) 

where 

{S?h, ),  =  RTbl  (Srbl  )h  +  SRTbl  (A  rbt  )h  (2.5-14) 

The  next  step  is  to  transform  the  body-object  displacement  vector  from  local 
level  coordinates  to  Earth  coordinates,  as 

(Arfc).«j£(ArJ,  (2.5-15) 

Errors  in  body  position  (latitude  and  longitude)  cause  errors  in  this  transforma¬ 
tion  and  thus  add  errors  to  the  body-object  displacement  vector  in  Earth  coor¬ 
dinates. 


(AiL),=(A rbl)e+(Srbl)e 


(2.5-16) 


where 

(K).  (2.5-17) 

The  last  step  is  to  add  the  body  position  vector  to  the  body-object  displacement 
vector  to  obtain  the  object  position  vector,  as 

(rt)e=(rh)e+(Arb,)e  (2.5-18) 

Errors  in  body  position  (including  altitude)  also  contribute  at  this  step,  and 

(2.5-19) 

where 

{SfX^iSfXHSf^  (2.5-20) 


2.6  Definitions 

All  these  items  are  defined  above,  but  they  are  listed  here  for  reference.  A  few 
items  are  in  more  than  one  category. 

2.6.1  Points 

O  is  the  center  of  the  Earth  (the  center  of  the  Earth  model  ellipsoid). 

B  is  center  of  the  body  frame,  also  called  the  center  of  the  ship’s  INS. 
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S  is  the  center  of  the  sensor  system. 

T  is  the  center  of  the  aerospace  object,  the  point  on  the  object  sensed  by  the 
sensor. 

The  “centers"  of  the  body  frame,  the  sensor  system,  and  the  object  are  not 
geometric  centers,  but  merely  points  fixed  to  the  items  in  question. 

2.6.2  Displacement  vectors 

rb  is  the  body  coordinate  frame  position  vector  that  goes  from  O  to  B . 

7,  is  the  aerospace  object  position  vector  that  goes  from  O  to  7  . 

A rbs  is  the  body-sensor  lever  arm  vector  that  goes  from  B  to  S  . 

A rst  is  the  sensor-object  displacement  vector  that  goes  from  S  to  T . 

A rbt  is  the  body-object  displacement  vector  that  goes  from  B  to  T, 

2.6.3  Coordinate  systems 

e-frame,  Earth  coordinates 

1-frame,  local  level  coordinates  (at  the  body  coordinate  frame  position) 
b-frame,  body  coordinates,  fixed  to  the  ship 
s-frame,  sensor  coordinates,  fixed  to  the  sensor 

2.6.4  Coordinate  transformation  matrices 

Rle  transforms  from  Earth  to  local  level  (at  the  INS)  coordinates. 

Rbl  transforms  from  local  level  (at  the  INS)  to  body  coordinates. 

Rbs  transforms  from  sensor  to  body  coordinates. 

2.6.5  Measured  items 

(0,A,h)  are  the  geodetic  coordinates  of  the  INS  ,  which  are  the  geodetic  latitude, 
the  longitude,  and  the  altitude  of  the  INS  above  the  ellipsoid. 

( y/h,0b,<pb )  are  the  three  Euler  angles  that  describe  the  body  coordinate  frame 

orientation  with  respect  to  local  level.  They  are  the  body  heading,  the  body 
pitch,  and  the  body  roll. 

(A rhs)b  is  the  body-sensor  lever  arm  vector  in  body  coordinates. 

Rbs  is  the  transformation  matrix  from  sensor  to  body  coordinates. 
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(p,y/s,6s)  are  the  sensor  measurements  of  range,  azimuth  {or  bearing),  and  ele¬ 
vation. 

In  all,  there  are  fifteen  scalar  error  sources  associated  with  the  fifteen  meas¬ 
ured  quantities.  Of  course,  beyond  these  fundamental  measured  items,  there 
are  many  derived  items.  Associated  with  each  of  the  fifteen  measured  quanti¬ 
ties  is  an  error.  These  errors  are: 

•  Latitude 

•  Longitude 

•  Height 

•  Three  navigation  tilts 

•  Three  lever  arm  errors 

•  Sensor  range 

•  Sensor  bearing 

•  Sensor  elevation 

•  Three  sensor  alignments 
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2.7  The  Aerospace  Object  Vector  Processing  Chain 

While  the  entire  chain  of  aerospace  object  vector  processing  has  been  dis¬ 
cussed  above,  with  error  considerations,  it  may  be  useful  to  briefly  summarize 
that  chain  from  a  strictly  nominal  (error  free)  point  of  view.  All  the  symbols 
here  are  defined  above,  and  all  the  equations  have  been  previously  presented. 
The  starting  point  is  the  sensor  measurements  of  the  object,  and  the  ending 
point  is  the  object  position  vector,  which  is  the  vector  from  the  center  of  the 
Earth  to  the  aerospace  object,  expressed  in  Earth  coordinates,  Figure  2-6. 


Step  1 :  Sensor  system  measures  the 
sensor-object  displacement  vector  in 
sensor  coordinates 


The  object  position  vector  in 
WGS  84  coordinates 


Step  2:  Sensor-object  displacement 
vector  is  transformed  to  body 
coordinates 


Step  3:  Body-sensor  level  arm  is  added 
to  produce  the  body-object  displacement 
vector  in  body  coordinates 


Step  4:  The  body-object  displacement 
vector  is  transformed  to  local  level 
coordinates 

Step  5:  The  body-object  displacement 
vector  is  transformed  to  Earth  coordinates 


Step  6:  The  navigation  system  position  vector 
is  added  to  produce  the  object  position  vector 
in  WGS  84  coordinates 


Figure  2-6  The  Aerospace  Object  Vector  Processing  Chain 

Step  1:  The  sensor  system  measures  the  sensor-object  displacement  vec¬ 
tor  in  sensor  coordinates 


(&*),  =P 


Ct 9, 


(2.7-1) 


Step  2:  The  sensor-object  displacement  vector  is  transformed  to  body  co¬ 
ordinates 
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(ArJb  =  RhMrsl)s 


(2.7-2) 


Step  3:  The  body-sensor  lever  arm  is  added  to  produce  the  body-object 
displacement  vector  in  body  coordinates 

(Arht  )b  =  (A  rbs  )h  +  (A  r„  )b  (2.7-3) 

Step  4:  The  body-object  displacement  vector  is  transformed  to  local  level 
coordinates 


(A  rb!)i  —  Rh<  (Arb,  )b 


(2.7-4) 


Step  5:  The  body-object  displacement  vector  is  transformed  to  Earth  co¬ 
ordinates 


(Art,  ),=/£(  Arfc), 


(2.7-5) 


Step  6:  The  navigation  system  position  vector  is  added  to  produce  the 
object  position  vector  in  WGS  84  coordinates 

(r,)e=(fb)e  +  (Arte)e  (2.7-6) 

These  same  steps  are  used  to  derive  measured  items  from  the  fundamen¬ 
tal  error  free  items. 

The  chain  of  errors  is  provided  by  the  following  summary  of  previously 
developed  formulas. 

Step  1:  Error  in  the  sensor-object  displacement  vector  in  sensor  coordi¬ 
nates 

(AFJ,  =(Arst)s  +(<&<, )*  (2.7-7) 

Cgs  C¥i  $8,  Cys  ~  C 6,  5V, 

{Srs,)s=Sp  C,s¥t  ~pdesSeS¥)  +p6yrs  C0CV>  (2.7-8) 

-Sg  Ce  0 

L  e>  Ji  L  Ji  L  Js 

Step  2:  Error  in  the  sensor-object  displacement  vector  in  body  coordi¬ 


nates 


<Ara)t  =(Ar„),+(«:„)( 
(ff„)s=iU<5r„  ),+*„(  Ar„), 


bs  bs 


(2.7-9) 

(2.7-10) 

(2.7-11) 
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(2.7-12) 


0 

~5y/2~ 

~<V," 

-Sy/, 

0 

Sy/{ 

—  SD.h 

<=>  (<¥)*  = 

Sy/2 

Sy/2 

-fy\ 

0 

Step  3:  Error  in  the  body-object  displacement  vector  in  body  coordinates 

(AFftf),=(A rb, )„+(&„,)„  (2.7-13) 

(K)b=(K)„  +  (<*A  (2.7-14) 

Step  4:  Error  in  the  body-object  displacement  vector  in  local  level  coor¬ 
dinates 


(Arhl),  =  (Arbt)l+(Srbt), 
(Srbl)l=RTbl(Srbl)b+SRlt(Arbl)b 


SRhl  -  &Qb,Rbl 


0 

-seD 


soL 

o 


-seE 

ddN 


(2.7-15) 

(2.7-16) 

(2.7-17) 

(2.7-18) 


mb  = 


l 

se,  - 

5dN 

o  J 

T 

0 

II 

£ 

0 

+sob 

+S¥b 

Ce>S<h 

b 

0 

b 

_~V 

b 

CeC^ 

(2.7-19) 


Step  5:  Error  in  the  body-object  displacement  vector  in  Earth  coordi¬ 


nates 


(Arbt)e^(Arbl)e+(SFbl)e 

(Kl-RU&J'+tel  (Arto), 


&,^foleRle 


<*4  = 


0 


3 

0 


&j>2  -5<px 


-S02 

0 


(2.7-20) 

(2.7-21) 

(2.7-22) 

(2.7-23) 


2-20 


7.2.6.3_TR(4-0Q2)_1.0_JSSEO_040407 


Step  6:  Error  in  the  aerospace  object  position  vector  in  Earth  coordinates 


(*))«. 


(<Sr*),  =  (rm  +  ^)<50 

-s,cx' 

-V* 

+  (rew  +  h)CfSA 

... . J 

Sh 

‘c#c/ 

e 

L  0  J 

e 

.  . 

(2.7-25) 

(2.7-26) 

(2.7-27) 
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2.8  Time  Error 


In  the  development  thus  far,  time  has  not  entered  the  problem.  The  point  of 
view  that  was  taken  is  that  there  is  a  single,  instantaneous,  measurement  of 
the  aerospace  object  by  the  sensor,  accompanied  by  supporting  information — 
at  that  same  instant  of  time — about  the  INS  position  and  attitude,  the  body¬ 
sensor  lever  arm,  and  the  relationship  of  the  body  coordinate  frame  to  the  sen¬ 
sor  coordinate  frame.  That  single  measurement  set  is  then  converted  into  a 
measured  object  position  with  respect  to  the  Earth.  However,  in  reality  there 
are  peer-to-peer  and  peer-to-reference  time  errors. 

A  better  conceptual  model  for  what  typically  happens  is  that  a  sequence  of 
such  sensor  measurements  (each  accompanied  by  the  time-dependent  support¬ 
ing  information)  is  made.  Each  of  these  measurements  results  in  a  computed 
object  position  in  Earth  coordinates.  The  sequence  of  object  positions  is  then  fit 
to  some  object  trajectory  model,  such  as  a  constant  velocity  model  or  a  con¬ 
stant  acceleration  model,  as  appropriate,  to  estimate  the  trajectory  as  a  func¬ 
tion  of  time.  Instead  of  resulting  in  just  a  single  object  position,  the  result  is  an 
object  position  vector  at  some  specified  time,  along  with  a  object  velocity  vector 
at  that  same  time,  and  possibly  a  object  acceleration  vector  at  that  time,  as 
well  (see  Figure  2-7).  These  three  vectors  are  accompanied  by  a  sensor  time  tag 
indicating  the  time  at  which  they  are  valid. 

The  errors  in  any  one  estimate  of  object  position  have  been  discussed  in  the 
previous  sections.  The  way  those  errors  combine  to  produce  errors  in  the  tra¬ 
jectory  estimate  is  not  considered  here.  That  process  depends  on  the  details  of 
the  trajectory  model  used,  the  fitting  process  used,  and  especially  on  the  time 
correlation  properties  of  the  errors  in  the  individual  measurements  and  their 
supporting  data. 

The  object  trajectory  information  is  then  usually  reported  to  some  other  par¬ 
ticipant,  possibly  remote,  to  use  in  dealing  with  the  object.  Let  us  call  that 
other  participant  the  user.  Typically,  the  user  is  not  interested  in  the  object  po¬ 
sition  at  the  time  it  was  sensed  or  reported,  but  is  interested  in  the  object  posi¬ 
tion  (and  perhaps  velocity)  at  some  later,  extrapolated,  time  called  the  time  of 
use.  The  object  position  derived  by  the  user  has  errors  for  three  reasons. 

1.  Even  if  all  the  information  about  the  object  trajectoiy  from  the  sensor  is  ab¬ 
solutely  correct,  there  will  still  be  errors  in  extrapolating  that  information  to 
the  time  of  use.  These  errors  are  called  extrapolation  errors. 

2.  The  object  position,  velocity,  and  acceleration  data  provided  by  the  sensor  is 
in  error,  and  these  errors  will  cause  errors  in  the  user-extrapolated  items  at 
a  later  time.  These  errors,  both  the  original  data  errors  of  the  sensor  and  the 
resulting  errors  by  the  user,  will  be  referred  to  as  trajectory  data  errors. 
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3.  The  time  tag  provided  by  the  sensor  is  in  error  (with  respect  to  some  ideal¬ 
ized  absolute  time  scale),  and  the  user  clock  also  has  errors.  The  difference 
of  these  errors  is  a  time  error,  which  will  cause  errors  in  the  extrapolated 
trajectoiy  items  by  the  user.  These  errors  in  the  user  items  will  be  called 
time  synchronization  errors  here.  All  three  types  of  errors  combine  to  pro¬ 
duce  what  are  lumped  here  as  “timing  errors”. 


t  —  Time  aerospace  object  is  observed  by  the  radar 

ro  =  Object  position  at  to 
vo  =  Object  velocity  at  time  to 
ao  =  Object  acceleration  at  time  to 

t  -  Time  object  track  data  is  used 

At  =  Nominal  time  latency  or  nominal 
time  past  to  target  track  data  is  used 

dt  =  t-(t0  +  A t)=  Time  error 

r(t0  +  At)  ~  Object  position  at  error-free  time  to  +At 
dr=~rit0+At  +  <5t)-r(t0  +  At)=  Track  error 


Figure  2-7  Track  Error  Due  To  Time  Error 


Trying  to  consider  all  three  types  of  errors  at  once  can  be  confusing,  since  they 
interact.  Additional  confusion  occurs  because  time  is  represented  as  the  inde¬ 
pendent  variable  in  the  error  equations,  so  time  errors  enter  in  a  different  way 
than  do  the  other  errors.  The  goal  here  is  a  linearized,  first-order  error  analysis, 
meaning  that  the  results  are  only  valid  to  first  order  in  the  (presumed  small) 
error  quantities.  Interactions  among  various  errors  may  be  ignored,  and  each 
error  may  be  treated  as  if  the  others  did  not  exist,  and  the  results  of  the  three 
treatments  simply  summed. 

The  vector  notation  of  the  preceding  development  is  appropriate  for  the  situa¬ 
tion  when  several  coordinate  systems  are  involved  and  everything  happens  at 
the  same  time.  It  is  not  so  appropriate  for  this  situation,  where  all  vectors  are 
expressed  in  one  coordinate  system — Earth  coordinates— and  time  is  an  issue. 
Here,  it  is  important  to  denote  the  time  dependence  and  not  important  to  de¬ 
note  the  coordinate  system  dependence.  The  focus  here  is  always  on  the  aero- 
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space  object,  so  it  is  not  necessary  to  subscript  trajectory  items  to  indicate  the 
object.  The  notation  will  change  to  reflect  these  differences  in  analysis. 

The  object  position  (and  velocity  and  acceleration)  are  reported  in  an  absolute 
positional  reference  system.  The  same  is  true  of  time.  That  is,  there  is  some  ab¬ 
solute  time  scale — such  as  universal  coordinated  time  (UTC),  or  GPS  time — 
used  to  report  time.  With  respect  to  this  scale,  both  the  sensor  and  the  user 
clocks  are  in  error. 

2.8.1  Extrapolation  errors 

The  (true)  time  at  which  the  sensor  reports  of  the  trajectory  items  are  valid  will 
be  denoted  by  t0 ,  known  as  time  zero.  The  true  position,  the  true  Earth-relative 

velocity,  and  the  true  Earth-relative  acceleration  of  the  aerospace  object  (all  ex¬ 
pressed  in  Earth  coordinates),  as  functions  of  the  (true)  time  will  be  denoted  as 
r(t) ,  v(r) ,  and  a(t ) .  For  convenience,  the  values  of  these  vectors  at  time  zero 

will  be  denoted  as  r0 ,  v0 ,  and  a0 . 

The  basic  relationships  among  these  quantities  are  that  velocity  is  the  time  de¬ 
rivative  of  position,  and  acceleration  is  the  time  derivative  of  velocity.  These  re¬ 
lationships  are  more  usefully  expressed  in  integral  terms  as 


v(t)  =  v0  +  J  a{T)dt 

(2.8-1) 

r(0  =  r0  +  |  v(z)dt 

(2.8-2) 

r(t)~  r0  +v0(r-i0)+  |  fd(A)dAdr 

(2.8-3) 

If  the  object  acceleration  were  constant  (in  general  it  is  not),  then  it  would  be 
true  that 

v(t)  =  v0+aQ(t-~t0)  (2.8-4) 

r{t)  =  r0  +v0(t-t0)  +  ~a0(t~t0)2  (2.8-5) 

Since  equations  of  this  form  are  assumed  to  be  used  to  extrapolate  the  object 
velocity  and  position,  let  us  define  the  “constant  acceleration”  items 

va(0  =  v0+50(r-/0)  (2.8-6) 

ra  (0  =  r0  +  v0  (t  - 10 )  +  ^  a0  (t  - 10 ) 2  (2.8-7) 


Then  we  may  write 
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(2.8-8) 


v(0  =  v,(0  +  J[5(r)-a0]rfr 

r(t)  =  ra (/)  +  f  f[a(yl)~a0]aUJr  (2.8-9) 

and  denote  the  differences  between  the  actual  items  and  the  constant  accelera¬ 
tion  items  by 

<^(0  =  v(0-va(0=  |[a(r)~a0]</r  (2.8-10) 

Srx(t)  =  r(t)-ra(t)=  j^[a(Z)~a0]dJUT  (2.8-11) 

If  all  the  trajectory  information  provided  by  the  sensor,  including  the  time  tag, 
were  absolutely  correct,  and  if  the  user  clock  were  also  error  free,  the  values  of 
the  two  differences  above  (evaluated  at  the  time  of  use)  would  be  the  errors  in 
the  extrapolated  velocity  and  position.  Thus  the  extrapolation  errors  are 

8*x  =  S5X (L )  =  v(*„ )  -  va (tu )  -  \\a{T)  -  a0  ]rfr  (2.8- 1 2) 

SFX  =  drx  (i tu )  =  r{tu  )~ra(tu)=  f‘J[aU)-a0]dAdT  (2.8-1 3) 

2.8.2  Trajectory  data  errors 

The  three  vector  trajectory  items  provided  by  the  sensor,  with  errors,  are  (defin¬ 
ing  new  notation) 

r0  =  r0  +  Sr0  =  r  (/0 )  =  r(t0 )  +  dr(t0 )  (2.8- 14) 

v0  =  i?0+<£0  =  v  (r0 )  ==  v(r0 )  -h  <S7(f0 )  (2.8-15) 

a0  -a0+Sa0  =  a(t0)  =  a(t0)  + &(tQ)  (2.8-16) 

Deviations  of  the  object  trajeetoiy  from  the  constant  acceleration  model  cause 
the  extrapolation  errors.  In  treating  the  trajectory  data  errors  we  assume  that 
the  constant  acceleration  model  is  correct.  The  error  is  described  by  equations 
2.8-14  to  2.8-16. 

Let  the  extrapolation  time  difference  or  network  latency  be  denoted  as 
A t  =  tu-t0.  Using  equations  2.8-1  to  2.8-3  with  constant  acceleration  leads  to 

the  resultant  trajectory  data  errors  given  by 

&d  =  va(tu)~va(tu)  =  &0+&0At  (2.8-17) 
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(2.8-18) 


Srd  =  K  (tu  )-ra(tu)  =  SrQ  +  dv0At  +  ^&0At2 

2.8.3  Time  synchronization  errors 

The  time  tag  reported  by  the  sensor  will  be  denoted  as  t0 ,  where  t0~t0+  <5r  v , 
with  dts  being  the  time  error  of  the  sensor,  see  Figure  2-8.  The  time  of  use,  on 
the  user’s  clock,  is  measured  as  7u ,  where  7U  =  tu  +  dtu ,  and  Stu  is  the  user  time 
error. 

The  focus  here  is  on  the  errors  made  in  extrapolating  the  trajectoiy  estimate  to 
a  position  estimate  at  some  future  time  point,  and  on  the  errors  caused  by  time 
synchronization  errors  between  the  sensor  clock  and  the  user  clock.  There  may 
be  other  time  synchronization  errors.  The  vehicle  navigation  system  might  have 
its  own  clock  that  differs  from  the  sensor  system  clock,  even  though  both 
clocks  are  nearby  on  the  same  vehicle.  These  other  possible  timing  errors  are 
not  considered  here. 

The  “measured”  time  for  extrapolation  is 

At  =tu -70=tu  +  &u -t0-Sts  =  A  t  +  St  (2.8-19) 

where  the  net  time  synchronization  error  is  St  —  ftu  —  <S5 .  That  is.  At  is  the  time 

actually  used  by  the  user  to  extrapolate  the  sensor  trajectoiy  data,  instead  of 
using  the  true  network  latency  interval  At .  The  errors  in  the  extrapolated  tra¬ 
jectory  items  produced  by  the  time  synchronization  error  are 

<5?,  =  a0St  (2.8-20) 

8rt  =  v0<£  +  50A/<5r  (2.8-21) 


2-26 


7.2.6.3_TR(4-002)_i  ,0_JSSEO_040407 


bt  =  Tu-Ts  =  btu-bts 


-4 - j. . 

btu  —  Tu  t 

user  i  f 

|  t . 

■  I  1  T 

measured  lime 

!  1 

1  1  lu 

-2  -1  0  j  1 

2  3 

sensor  ,  , 

i 

1  1  |  r 

measured  time  '  ' 

-2  -I 

0 

1  2  3 

true  time 

] 

I  |  . | . 

(UTC(USND))  1  1 

1  1  I 

-2  -1  0  12  3 


Figure  2-8  Time  Synchronization  Error 

2.8.4  Combined  errors 

Putting  the  three  results  together— by  simply  adding  them — yields  the  com¬ 
bined  error  in  the  trajectoiy  estimates  by  the  user  as 

8Se  =  5$x  +  +&,  =  |“[a(r)  -  a0  \lx  +  <*0  +  <5a0A/  +  20<5r  (2.8-22) 

<3rc  =  5rx  +  8rd  +  Srt  (2.8-23) 

Src  =  f  f  [5(A)  -  50 +  <5r0  +  #0Af  +  l&0At2  +  v0£  +  50Ar£  (2.8-24) 

"0  "o  2 

Recapping,  we  have  the  following  interpretation  for  the  terms  on  the  right  side 
of  equation  2.8-24: 

•  I"  J *\a{X)-aa\lMT  is  the  error  due  to  using  an  average  acceleration 
in  place  of  the  true  acceleration 

•  Sr0  is  the  error  in  the  aerospace  object  position  vector  in  Earth  co¬ 
ordinates  given  by  equation  2.7-26 
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•  Sv0At  +  -Sa0At2  represents  the  error  in  the  velocity  and  acceleration, 

respectively,  used  in  the  tracking  algorithm  to  extrapolate  the  ob¬ 
ject’s  position 

•  v0$  +  a 0  At  St  is  the  contribution  to  track  error  due  to  the  total  time 
error. 

2.8.5  The  processing  sequence 

At  the  risk  of  belaboring  the  issue,  the  sequence  of  processing  steps  by  the 
user  will  be  outlined.  The  input  data  to  this  sequence  are  the  trajectory  items 
and  the  aceompaning  time  tag  provided  by  the  sensor.  That  is,  the  inputs  are 
(F0,v0,50,ro).  One  additional  input,  provided  by  the  user,  is  the  measured  time  at 
which  the  trajectory  estimate  is  to  be  used,  or  tu .  The  user  computes  the 
measured  extrapolation  time  difference  as  A 7  =  tu  —t0,  and  then  extrapolates  the 
trajectory  items  using  the  constant  acceleration  model 


v{tu)  =  v0+a0At 

~  _  ~  ~  ~  \~  ~2 
r(tu)  =  r0+v0At  +~a0At 

The  errors  in  this  process  are  given  by 

v(tu)~v(tu)  =  &c 

r(tu)-r(tu)  =  Src 

where  the  combined  errors  are  delineated  above. 


(2.8-25) 

(2.8-26) 

(2.8-27) 

(2.8-29) 
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2.9 


The  Sensitivities 


The  analysis  presented  in  the  previous  sections  is  characterized  as  a  linear  er¬ 
ror  analysis  since  only  error  terms  up  to  first  order  are  retained.  As  such  the 

results  can  be  expressed  in  the  compact  form  [$rt\=Se ,  where  S  is  the  sensitiv¬ 
ity  matrix  and  e  is  the  error  vector  defined  by 


North  Position  Error 
East  Position  Error 
Height  Error 
Navigation  Tilt ,  N 
NavigationTilt,  E 
Navigation  Tilt,  D 
Lever  Arm  Error,  £>, 
Lever  Arm  Error,  h2 
Lever  Arm  Error,  b. 
Sensor  Range  Error 
Sensor  Bearing  Error 
Sensor  Elevation  Error 
Sensor  Misalignment,  yx 
Sensor  Misalignment,  lf/2 
Sensor  Misalignment,  l//2 
Time  Error 


(2.9-1) 


This  form  is  advantageous  since  the  covariance  of  track  error  can  be  directly 
obtained  from  the  expression  cov(Sr)~Scov(e)ST .  Similarly,  the  covariance  of 
the  track  error  for  groups  of  error  contributors  can  be  calculated  by  consider¬ 
ing  the  appropriate  submatrices  of  S. 

The  elements  of  S  provide  information  on  how  fast  the  components  of  track  er¬ 
ror  change  as  the  individual  error  terms  change.  The  three  rows  of  S  corre¬ 
spond  to  the  x,  y,  and  z  components  of  the  WGS  84  coordinate  system  and 
each  column  of  S  corresponds  to  an  element  of  the  error  vector.  For  example, 
we  will  denote  S(*,  Sensor  Range  Error)  as  the  three  sensitivities  (x,  y,  and  z  di¬ 
rections,  respectively)  corresponding  to  the  sensor  range  error. 


The  elements  of  the  sensitivity  matrix  are  obtained  by  rearranging  the  expres¬ 
sions  provided  in  Section  2.6  and  making  the  following  two  identifications. 
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First,  latitude  error  is  equivalent  to  north  position  error  and  second,  east  posi¬ 
tion  error  is  obtained  by  dividing  longitude  sensitivity  by  the  cosine  of  latitude. 
In  what  follows,  expressions  of  the  form  A ra/t  represent  the  1th/  element  of  the 

displacement  vector  from  a,  where  a  represents  the  sensor  or  body,  to  the  ob¬ 
ject,  expressed  in  the  j8-frame.  The  columns  of  S  are  given  by  the  following: 


S(*,  North  Position  Error)  =  (rm+h) 


“Vi" 

I 

1 

> 

_ 1 

“Vi 

0 

- j 

O* 

_ \ 

e 

- 1 

< 

_ i 

(2.9-2) 


S(*,East  Position  Error)  =  (rn  +h) 


-V 

S*Ar»‘n 

cx 

+  (Rl/Cfi) 

~  V*,,  ~C*Arb,n 

0 

e 

C#Arbrn 

(2.9-3) 


S(*,  Height  Error)  = 


C,c, 


(2.9-4) 


S(*,  Navigation  Tilt,  N)  -  RfeRl, 


0 

-Ar* 

Ar„. 


(2.9-5) 


S(* , Navigation  Tilt,  E)  -  RjeRh 


Ah< 

0 

-An, 


J b 


(2.9-6) 


S(*, Navigation  Tilt,  D)  =  RjeRh 


-Arb, 

Ar,„ 

0 


(2.9-7) 


T  dT 


S(*, Lever  Arm  Error ,  bx)  =  Rle R 


(2.9-8) 


r  »t 


S{*, Lever  Arm  Error,  b2 )  =  Rk,  R, 


(2.9-9) 
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S(*t Lever  Ann  Error,  b})  =  RkRt 


T  r>T 

bl 


S{*, Sensor  Range  Error)  =  RjeRl,Rbs 


C e,  Cy, 


S(*, Sensor  Bearing  Error)  =  RlRTMRbsP 


'  C(>,  S¥s 
cec¥t 
0 


S(*, Sensor  Elevation  Error)  =  RlRThlRbsp 


-S'C* 

-C 


S(*,  Sensor  Misalignment,  y/x )  =  RjeRT 


0 

Ar, 


«b) 


■A  r 


S(*,  Sensor  Misalignment,  y/2 )  =  Rle  R 


rDr 

bl 


Ar, 


0 

Ar, 


f  -  Ar, 


S(*,  Sensor  Misalignment, 


Ar. 


*»»i 


0 


S(*,Time  Error)  = 


Kov 


"0. 


+  A/ 


Ov 


*0, 


(2.9-10) 


(2.9-11) 


(2.9-12) 


(2.9-13) 


(2.9-14) 


(2.9-15) 


(2.9-16) 


(2.9-17) 
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3.  RESULTS 


3. 1  Overview 

It  should  come  as  no  surprise  that  the  analysis  results  depend  strongly  on  the 
particular  scenario,  the  operational  conditions,  and  the  nominal  error  budget. 
The  Common  Reference  Scenarios  (CRS)  are  being  developed  to  provide  a  uni¬ 
form  set  of  threat  and  operational  scenarios  for  analysis  purposes.  When  these 
CRS  are  completed  much  of  the  scenario-dependent  variability  will  be  elimi¬ 
nated.  However,  at  the  time  of  this  analysis  the  CRS  were  not  finalized  and  we 
have  chosen  parameter  values  that  represent  realistic  conditions  for  a  sea- 
based  system  in  order  to  identify  the  important  sensitivities  and  to  quantify  the 
relative  importance  of  the  error  contributors. 

We  will  use  the  numerical  values  of  the  sensitivity  matrix  to  illustrate  how  the 
various  errors  are  mapped  into  track  error  and  provide  insight  into  the  role  the 
geometry  of  the  problem  plays.  When  error  budgets  are  discussed  we  will 
measure  track  accuracy  in  terms  of  Spherical  Error  Probable  (SEP),  which 
represents  the  radius  of  a  sphere  centered  at  the  nominal  location  of  the  aero¬ 
space  object  that  encompasses  50%  of  the  error  distribution. 

3.2  Nominal  Scenario  and  Operational  Parameter  Values 

The  parameter  values  for  the  operational  scenario  are  provided  in  Table  3- 1 . 
These  parameters  place  the  ship  in  the  Atlantic  Ocean  sensing  an  object  low  in 
elevation  at  a  range  of  100  nm.  The  aerospace  object  represents  the  capabilities 
of  a  high-speed  aircraft.  The  roll,  pitch,  and  yaw  of  the  ship  are  set  to  zero  for 
convenience,  without  loss  of  generality. 


3-1 


7.2.6.3_TR(4-002)_1.0_JSSEO_040407 


Table  3-1  Nominal  Operational  Parameter  Values 


Item 

Parameter 

Units 

Value 

Latitude 

deg 

45 

Longitude 

deg 

-70 

Neiqht 

m 

0 

Roll 

deg 

0 

Ship 

Pitch 

deg 

0 

Yaw 

deg 

0 

Lever  arm, 

m 

-30 

Lever  arm,  b2 

m 

2 

Lever  arm,  b3 

m 

8 

Range  to  target 

nm 

100 

Sensor 

Bearing  to  target 

deg 

45 

Elevation  to  target 

deg 

2 

Network 

Latency 

sec 

5 

Velocity,  x 

kt 

500 

Velocity,  y 

kt 

0 

Aerospace 

Velocity,  z 

kt 

0 

Object 

Acceleration,  x 

9 

9 

Acceleration,  y 

g 

0 

Acceleration,  z 

_ 2 _ 

0 

3.3  Sensitivity  Matrix 

The  transpose  of  the  sensitivity  matrix  is  provided  in  Table  3-2.  The  numerical 
entries  in  the  table  represent  the  amount  of  track  error  in  each  of  the  three 
axes  of  the  WGS  84  coordinate  frame  that  are  generated  by  a  unit  error  term. 
For  example,  if  x  is  the  first  dimension  of  track  accuracy  and  east  position  is 
the  second  error,  then  a  value  of  0.91 1  for  the  (1,2)  element  of  the  sensitivity 
matrix  means  that  eveiy  meter  of  east  position  error  causes  0.91 1  m  of  track 
error  in  the  positive  direction  of  x.  If  the  east  position  error  is  1 ,000  meters 
then  its  contribution  to  track  error  is  91 1  m  in  the  positive  x  direction. 
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Table  3-2  Sensitivity  Matrix  (Transposed) 


Sensitivities  (m/unit) 

Item 

Error  Term 

Unit 

WGS-84 

Earth  Fixed  Axes 

X 

y 

z 

North  Position 

m 

-0.24706 

0.67879 

0.69329 

East  Position 

m 

0.91149 

0.36258 

0 

Height 

m 

0.24184 

-0.66446 

0.70711 

Navigation  tilt,  N 

mrad 

-25.58618 

89.17173 

-92.545 

Ship 

Navigation  tilt,  E 

mrad 

33.20572 

-91.23196 

87.95773 

Navigation  tilt,  D 

mrad 

154.60765 

-42.21175 

-92.545 

Lever  arm  error,  pi 

m 

-0.24184 

0.66446 

0.70711 

Lever  arm  error,  p2 

m 

0.93969 

0.34202 

0 

Lever  arm  error,  p3 

m 

-0.24184 

0.66446 

-0.70711 

Range  to  target 

m 

0.50159 

0.68807 

0.52437 

Bearing  to  target 

mrad 

154.63536 

-42.20016 

-92.54359 

Elevation  to  target 

mrad 

41.57299 

127.58352 

127.64471 

Sensor 

Sensor  misalignment,  xpi 

mrad 

25.57818 

-89.17314 

92.54359 

Sensor  misalignment,  ipi 

mrad 

-33.21491 

91.25721 

-87.97329 

Sensor  misalignment,  ipa 

mrad 

154.63536 

42.20016 

92.54359 

Network 

lime  error* 

sec 

699.15174 

0 

0 

Table  3-3  Unit  Radial  Track  Error 


Sensitivities  (m/unit) 

UnitRa- 

Item 

Error  Term 

Unit 

WGS-84  Earth  Fixed  Axes 

dial 

X 

y 

z 

lxacK 

Error 

North  Position 

m 

-0.247 

0.679 

0.693 

1.0 

East  Position 

m 

0.911 

0.363 

0.000 

1.0 

Height 

m 

0.242 

-0.664 

0.707 

1.0 

Navigation  tilt,  N 

mrad 

-26 

89 

-93 

131.0 

Ship 

Navigation  tilt,  E 

mrad 

33 

-91 

88 

131.0 

Navigation  tilt,  D 

mrad 

155 

-42 

-93 

185.1 

Lever  arm  error,  pi 

m 

-0.242 

0.664 

0.707 

1.0 

Lever  arm  error,  p2 

m 

0.940 

0.342 

0.000 

1.0 

Lever  arm  error,  p3 

m 

-0.242 

0.664 

-0.707 

1.0 

Range  to  target 

m 

0.502 

0.688 

0.524 

1.0 

Bearing  to  target 

mrad 

155 

-42 

-93 

185.1 

Elevation  to  target 

mrad 

42 

128 

128 

185.2 

Sensor 

Sensor  misalignment,  ipi 

mrad 

26 

-89 

93 

131.0 

Sensor  misalignment,  i f2 

mrad 

-33 

91 

-88 

131.0 

Sensor  misalignment,  m 

mrad 

155 

42 

93 

185.1 

Network 

lime  error* 

sec 

699 

0 

0 

699.2 

•Represents  the  sensitivity  of  any  time  error  such  as  time  synchronization,  timestamp,  etc. 
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Insight  into  how  individual  errors  contribute  to  track  error  can  be  obtained  by 
examining  the  root-sum-square  (RSS)  of  the  elements  of  the  sensitivity  matrix. 
These  values  are  provided  in  the  column  titled  “Unit  Radial  Track  Error  in  Ta¬ 
ble  3-3. 

The  values  of  unity  for  Unit  Radial  Track  Error  corresponding  to  the  three  rows 
of  navigation  position  error  imply  that  a  meter  of  position  error  in  north  or  east, 
or  height,  translates  one-for-one  into  a  meter  of  track  error.  This  mapping  is 
independent  of  operational  conditions  and  scenario,  e.g.,  aerospace  object  con¬ 
ditions.  The  same  holds  for  the  contribution  of  lever  arm  errors,  and  sensor 
range  error. 

If  the  RSS  of  the  three  navigation  tilt  Unit  Radial  Track  Error  values  is  com¬ 
pared  to  the  RSS  of  the  three  sensor  misalignment  Unit  Radial  Track  Error  val¬ 
ues,  one  finds  that  the  difference  is  insignificant.  This  implies  that  the  contri¬ 
butions  of  these  two  misalignment  contributors  are  nearly  identical.  The  reason 
for  this  is  that  the  distance  separating  the  navigation  center  from  the  sensor 
center,  31  m,  is  small  compared  to  the  distance  from  the  sensor  to  the  aero¬ 
space  object,  100  run.  The  magnitude  of  the  contribution  of  a  navigation  or 
sensor  misalignment  error  is  proportional  to  the  distance  from  the  sensor  to 
the  object. 

The  contribution  of  sensor  bearing  error  to  track  accuracy  is  proportional  to 
range  to  the  object  times  the  cosine  of  the  object’s  elevation.  The  sensor  eleva¬ 
tion  error  contribution  to  track  accuracy  is  proportional  to  range  to  the  object. 
The  numerical  values  in  Table  3-3  are  almost  equal  because  the  cosine  of  2  de¬ 
grees  elevation  is  close  to  unity. 

The  sensitivity  of  the  time  contribution  is  dependent  on  the  aerospace  object’s 
velocity  and  acceleration  and  the  network  latency.  If  any  of  these  operational  or 
scenario  conditions  increase  (decrease)  then  the  sensitivity  will  also  increase 
(decrease). 

3.4  Nominal  Error  Values 

The  nominal  error  budget  for  the  navigation,  sensor,  and  time  errors  is  pre¬ 
sented  in  Table  3-4.  The  errors  are  all  assumed  to  be  normally  distributed  with 
zero  mean.  For  the  present  analysis  all  errors  are  assumed  to  be  independent. 
However,  nothing  in  the  error  analysis  of  Section  2  requires  this  to  be  true.  The 
uncorrelated  model  was  chosen  for  convenience.  More  detailed  analysis  of  the 
individual  subsystems  is  required  to  develop  the  more  general  correlated  error 
model. 
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Table  3-4  Nominal  Error  Budget 


Item 

Error  Parameter 

Unit 

One-Sigma  Value 

North  position 

m 

300 

East  position 

m 

300 

Heiqht 

m 

150 

Navigation  tilt,  N 

mrad 

0.1 

Ship 

Navigation  tilt,  E 

mrad 

0.1 

Navigation  tilt,  D 

mrad 

0.1 

Lever  arm  error,  b , 

m 

0.1 

Lever  arm  error,  b2 

m 

0.1 

Lever  arm  error,  b3 

m 

0.1  . 

Range  to  object 

m 

300 

Bearing  to  object 

mrad 

1 

Elevation  to  object 

mrad 

1 

Sensor 

Sensor  misalignment,  4^ 

mrad 

1 

Sensor  misalignment,  ip2 

mrad 

1 

Sensor  misalignment,  tp3 

mrad 

1 

Network 

Time  error* 

sec 

0.5 

•Represents  the  total  effect  of  all  time  errors  such  as  time  synchronization,  timestamp,  etc. 

The  error  budget  for  each  system  also  provides  the  statistical  description 
needed  to  introduce  the  proper  level  of  bias  into  each  system  associated  with  a 
platform  specific  version  of  the  Integrated  Architecture  Behavior  Model  (IABM). 
A  random  value  can  be  drawn  from  the  distribution  of  each  error  each  time  the 
complete  IABM  is  executed  either  during  testing  or  during  runs  within  the 
Joint  Distributed  Engineering  Plant  (JDEP)  technical  framework. 

3.5  Nominal  Track  Accuracy  Error  Budget 

Combining  the  sensitivity  matrix.  Table  3-2,  with  the  nominal  error  budget. 
Table  3-4,  leads  to  the  nominal  track  accuracy  error  budget  displayed  in  Figure 
3-1.  For  each  of  the  error  terms  listed  on  the  left  of  the  bar  chart,  the  length  of 
the  bars  shows  the  contribution  to  track  error,  expressed  as  Spherical  Error 
Probable  (SEP),  if  the  error  under  consideration  is  the  only  error  in  the  system. 
For  example,  if  north  position  error  were  the  only  error,  then  SEP  for  the  sys¬ 
tem  would  be  about  250  m. 

The  bar  chart  also  illustrates  the  cumulative  effect  of  several  errors.  The  indi¬ 
vidual  yellow  bars  roll  up  to  the  blue  bars,  the  individual  blue  bars  roll  up  to 
the  red  bars  and  the  three  red  bars  roll  up  to  the  total  error  bar  color-coded 
purple.  For  example,  the  sensor  range,  sensor  bearing,  and  sensor  elevation  er¬ 
ror  contributors  combine  to  give  a  total  sensing  error  contribution  of  approxi¬ 
mately  350  m  SEP.  The  sensor  misalignment  contribution  and  the  sensor 
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measurement  contribution  combine  to  yield  an  approximate  contribution  of 
425  m  SEP  for  the  sensor  subsystem.  The  three  major  error  groups,  navigation, 
sensor,  and  time,  combine  to  yield  a  total  system  error  of  approximately  650  m 
SEP. 


Error 

Total  Navigation 
Total  Position 
North  Position 
East  Position 
Height 
Total  Tilt  Error 
Navigation  Tilt  N 
Navigation  Tilt  E 
Navigation  Tilt  D 
Total  Lever  Arm 
Lever  Arm  1 
Lever  Arm  2 
Lever  Arm  3  P 
Total  Sensor 
Total  Sensing 
Sensor  Range 
Sensor  Bearing  p 
Sensor  Elevation 
Total  Sensor  Align 
Sensor  Alignment  1 
Sensor  Alignment  2 1- 
Sensor  Alignment  3  p 
Time 
Total  Error 


Aerospace  Object 


400 
SEP  (m) 


Figure  3-*l  Nominal  Track  Accuracy  Error  Budget 

The  error  budget  also  provides  a  quick  way  to  rank  the  importance  of  the  vari¬ 
ous  error  contributors  to  track  accuracy.  For  this  case,  the  short  length  of  the 
bars  representing  the  navigation  tilt  and  lever  arm  error  contributions  implies 
the  insignificance  of  these  error  terms.  That  is,  resources  would  be  wasted  on 
reducing  these  error  terms  since  their  overall  contribution  to  the  total  system 
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accuracy  is  so  small.  Any  further  reduction  in  these  errors  would  cause  a 
minimal  reduction  on  total  system  error. 

The  magnitude  of  the  contributions  also  suggests  that  when  testing  the  LABM 
the  error  terms  with  small  contributions  should  be  included  first.  Otherwise 
the  larger  error  terms  will  dominate  and  it  will  be  difficult  to  determine  if  the 
small  terms  are  being  implemented  correctly. 

The  nominal  system  also  represents  a  balanced  system  with  respect  to  naviga¬ 
tion  and  total  sensor  error.  The  contributions  from  these  two  error  categories 
are  nearly  equal.  This  implies  that  either  would  be  a  good  candidate  to  consider 
for  making  improvements  to  track  accuracy.  However,  which  is  the  easiest  to 
correct  is  another  matter  to  consider.  The  time  error  would  rank  last  among  the 
three  major  error  groups  as  a  candidate  for  improving  the  total  system  error 
because  its  contribution  is  relatively  small. 

3.6  Sensitivity  of  Track  Accuracy  to  Individual  Errors 

This  section  presents  a  number  of  charts  displaying  the  sensitivity  of  track  ac¬ 
curacy  to  changes  in  the  nominal  one-sigma  error  parameters.  These  charts 
provide  important  information  on  how  track  accuracy  changes  if  the  nominal 
values  change.  For  comparison  purposes  the  nominal  values  presented  in  Fig¬ 
ure  3- 1  correspond  to  where  the  light  gray  vertical  line  intersects  the  sensitivity 
curves  in  Figures  3-2  to  3-12.  This  sensitivity  information  is  relevant  if  there  is 
some  uncertainty  in  the  confidence  of  the  one-sigma  values  or  if  there  is  a  need 
to  make  predictions  of  what  the  total  system  accuracy  will  be  if  a  change  is 
made  in  the  one-sigma  value.  This  section  is  divided  into  three  subsections, 
each  devoted  to  one  of  the  three  major  error  groups,  navigation,  sensor,  and 
time.  The  displays  for  the  sensitivities  of  track  accuracy  to  time  are  presented 
in  a  different  format  to  show  the  effect  of  time  latency  on  system  accuracy  and 
to  cover  two  ways  in  which  time  error  can  be  mapped  to  system  level  track  er¬ 
ror. 

3.6. 1  Sensitivity  to  navigation  errors 

Figures  3-2  to  3-5  display  the  sensitivity  of  track  accuracy  to  changes  in  the 
one-sigma  values  for  north  position  error,  east  position  error,  height  error,  and 
navigation  misalignments,  respectively.  In  all  the  figures,  the  curve  represent¬ 
ing  the  sensitivity  of  the  position  contribution 
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Track  Accuracy  SEP  (m)  Track  Accuracy  SEP  (m) 


North  Position  One-sigma  (m) 

Figure  3-2  Sensitivity  to  North  Position  Error 


Figure  3-3  Sensitivity  to  East  Position  Error 
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T rack  Accuracy  SEP  (m)  Track  Accuracy  SEP  (m) 


Height  One-sigma  (m) 

Figure  3-4  Sensitivity  to  Height  Error 


Sensitivity  to  Navigation  Tltt 


Figure  3-5  Sensitivity  to  Navigation  Tilt 
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coincides  with  the  curve  representing  the  contribution  of  navigation  because 
the  misalignment  contribution  is  so  small.  Figures  3-2  and  3-3  show  that  the 
sensitivities  of  total  system  track  error  to  north  position  and  east  position  error 
are  very  nearly  equal.  The  sensitivity  is  less  pronounced  for  the  height  error. 
Figure  3-4,  and  is  very  small  for  the  total  navigation  misalignment  or  tilt.  Fig¬ 
ure  3-5. 

3.6.2  Sensitivity  to  sensor  error 

The  sensitivities  of  track  accuracy  to  the  sensor  errors  are  displayed  in  Figures 
3-6  to  3-1 1.  Unlike  the  navigation  subsystem  where  the  position  error  group 
dominated  the  navigation  contribution,  both  the  sensing  errors  (range,  bearing 
and  elevation)  and  the  sensor  misalignments  about  each  axis  are  important 
contributors  to  track  error.  Of  all  the  sensor  errors,  the  sensor  range  error  has 
the  largest  influence  on  track  accuracy.  This  is  followed  by  the  sensitivities  due 
to  bearing  error,  1^3-axis  misalignment,  elevation  error,  and  the  other  two  sen¬ 
sor  misalignments.  Because  the  system  is  balanced  in  terms  of  the  contribu¬ 
tions  to  track  error  from  the  navigation  and  sensor  subsystems,  decreases  in 
error  from  the  nominal  value  have  a  small  effect  on  reducing  the  total  system 
error.  Substantial  reductions  in  track  error  can  only  be  made  by  reducing  both 
the  navigation  position  and  the  sensor  errors. 
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T rack  Accuracy  SEP  (m)  Track  Accuracy  SEP  (m) 


Sensor  Range  One-sigma  (m) 

Figure  3-6  Sensitivity  to  Sensor  Range  Error 


Sensor  Bearing  One-sigma  (mrad) 


Figure  3-7  Sensitivity  to  Sensor  Bearing  Error 
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Track  Accuracy  SEP  (m)  Track  Accuracy  SEP  (m) 


1 000 


Sensitivity  to  Sensor  Elevation  Error 

1'  . . . . ,~T~’-T 


Sensor  Elevation  One-sigma  (mrad) 

Figure  3-8  Sensitivity  to  Sensor  Elevation  Error 


Sensor  Misalignment  One-sigma  (mrad) 

Figure  3-9  Sensitivity  to  Sensor  Misalignment,  t^i-Axis 
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Track  Accuracy  SEP  (m)  Track  Accuracy  SEP  (m) 


Sensor  Misalignment  One-sigma  (mrad) 

Figure  3-10  Sensitivity  to  Sensor  Misalignment,  t^-Axis 

SensHMty  to  Sensor  Misalignment  Prtfrtods 
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Figure  3-11  Sensitivity  to  Sensor  Misalignment,  1/^3-Axis 
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3.6.3  Sensitivity  to  time  error 

The  sensitivity  of  track  accuracy  to  time  error  is  displayed  in  Figure  3-12.  The 
curve  shows  that  reducing  the  time  error  to  zero  reduces  track  error  by  only  a 
little  more  than  50  m  SEP.  This  is  because  the  navigation  errors  and  sensor  er 
rors  dominate  the  error  budget.  Doubling  the  time  error  budget  so  it  is  bal¬ 
anced  with  navigation  and  sensor  will  increase  total  track  SEP  by  approxi¬ 
mately  100  m. 


Additional  insight  into  the  sensitivity  of  track  accuracy  to  time  error  is  obtained 
by  examining  the  sensitivities  for  an  expanded  set  of  threat  objects  and  net- 
work  latency  times.  The  velocity  and  acceleration  for  a  set  of  aerospace  objects 
is  provided  in  Table  3-5.  These  parameter  values  represent  maximum  capabili¬ 
ties  and  provide  worst  case  sensitivities  of  track  accuracy  error  to  time  error. 
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Table  3-5 


Aerospace 

Object 

Velocity 
Magnitude  (kt) 

Acceleration 
Magnitude  (g) 

Helicopter 

200 

3.5 

Aircraft  1 

300 

5 

Cruise  Missile 

500 

0.5 

Aircraft  2 

800 

9 

Aircraft  3 

1000 

9 

Aircraft  4 

1200 

9 

Aircraft  5 

1600 

9 

We  will  first  consider  how  the  time  errors  are  mapped  into  track  error.  As  dis¬ 
cussed  in  section  2,  the  sensitivity  of  track  error  to  time  error  is  given  by  the 

expression  vo  +  a0At  where  vo  is  the  object  velocity  at  the  time  the  object  is 

sensed,  t0 ,  and  ao  is  the  average  acceleration  over  the  nominal  latency  time, 

At  of  the  network.  We  will  first  consider  the  case  of  the  directions  of  velocity 
and  acceleration  for  the  object  as  being  parallel  and  show  the  relationship  be¬ 
tween  the  normal  densities  for  the  time  and  track  errors.  Figures  3-13  to  3-15 
illustrate  this  by  showing  the  mapping  of  the  one-sigma  values  for  time  error 
into  the  one-sigma  values  for  track  error  for  the  case  of  Aircraft  3  and  a  net¬ 
work  latency  of  5  sec  for  three  ranges  of  time  error.  For  the  case  of  Aircraft  3 
and  a  network  latency  of  5  sec,  the  sensitivity  is  956  m  of  track  error  for  eveiy 
sec  of  time  error.  We  see  that  time  errors  in  the  range  of  0.001  to  0.01  sec  lead 
to  track  errors  in  the  range  of  1  to  10  m,  time  errors  between  0.01  and  0.1  sec 
give  rise  to  track  errors  between  10  to  96  m,  and  for  the  third  case  of  time  er¬ 
rors  between  0.1  and  1  sec,  the  track  errors  are  between  96  and  956  m. 
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Velocity  magnitude = 1 ,000  kt  Acceleration  magnitude = 9  g 


We  will  next  show  the  sensitivities  of  track  accuracy  to  time  errors  for  each  of 
the  aerospace  objects  listed  in  Table  3-5.  So  far  we  have  assumed  that  the  ob 
ject  velocity  and  acceleration  are  known  and  that  these  vector  quantities  are  in 
the  same  direction.  Since  the  time  error  is  one-dimensional  and  the  directions 
of  the  velocity  and  acceleration  are  assumed  to  be  known,  the  density  of  track 
error  will  be  a  degenerate  trivariate  normal.  That  is,  the  density  will  be  one  di¬ 
mensional  along  the  direction  of  the  vector  quantities.  We  will  relax  the  as¬ 
sumption  of  known  directions  for  the  velocity  and  acceleration  by  considering  a 
second  case  where  we  let  the  directions  for  velocity  and  acceleration  be  inde¬ 
pendent  and  uniformly  distributed  in  three-space.  This  assumption  will  give 
rise  to  a  truly  three-dimensional  track  error  and  more  realistically  represent  a 
true  warfighting  scenario.  It  is  noted  that  we  will  still  consider  the  magnitudes 
of  the  velocity  and  acceleration  of  the  object  to  be  known.  Experts  in  the  field  of 
tracking  have  noted  that  prior  knowledge  of  the  object’s  velocity  and  accelera¬ 
tion  may  not  be  known  exactly  and  this  lack  of  knowledge  often  leads  to  track¬ 
ing  errors.  We  are  presently  investigating  how  to  model  this  situation  and  in¬ 
clude  it  as  a  tracking  error  in  future  error  budget  discussions. 

To  complete  the  present  discussion  we  will  examine  the  following  two  cases  pa¬ 
rametrically. 

•  Case  1 :  Assume  direction  and  velocity  and  acceleration  is  fixed  and 
known 

•  Case  2:  Assume  random  direction  of  velocity  and  acceleration  but 
magnitudes  are  known. 

We  will  parametrically  vary 

•  Magnitude  of  velocity  and  acceleration  by  considering  the  aerospace 
objects  presented  in  Table  3-5 

•  Nominal  network  time  latency.  At 

•  One-sigma  value  of  the  time  error. 

Track  accuracy  will  be  measured  by 

•  Standard  deviation  (one-sigma)  representing  69%  coverage  for  Case  1 

•  Spherical  Error  Probable  (SEP)  representing  50%  coverage  for  Case  2. 

The  full  set  of  sensitivity  curves  are  presented  in  Figures  3-16  to  3-23. 

Figures  3-16  and  3-17  show  the  sensitivities  of  track  accuracy  for  the  same 
helicopter  but  for  different  ranges  of  time  error.  Figure  3-16  provides  a  clear 
indication  of  the  magnitude  of  the  track  accuracy  ranging  from  a  minimum  of 
approximately  20  m  one  sigma  or  SEP  to  approximately  1,150  m  one-sigma 
and  825  m  SEP  for  time  errors  between  0.2  and  1  sec  and  for  all  cases  of  net- 
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work  time  latency.  Figure  3-17  narrows  in  on  the  0.02  to  0  1  sec  mnge ^or  time 
error  and  shows  the  largest  track  accuracy  Is  approximately  115  m  one-sigma 
for  Case  1  and  80  m  SEP  for  Case  2. 

Based  on  the  analytic  representation  of  the  sensitivity  and  the  family  of  sens!- 
tivity  curves  we  immediately  see  that  track  error 

•  Increases  as  the  velocity  of  the  object  increases 

.  *  Increases  as  the  acceleration  of  the  object  increases 

•  Increases  as  the  network  time  latency  increases. 

r„mnarintf  Figure  3-19  with  any  of  the  others  shows  the  combined  effect  of  ac- 

celeration^ancf  network  latency  on  track  accuracy^e  n^row  spread^bet^ 

hrcre^iesr^hen'kirger^accelerahor^are  combined  with  increased  network  time 

latency  -  a  multiplicative  effect  is  taking  place. 

F-iffure  3-23  shows  that  under  worst  case  conditions,  the  track  accuracy  can  be 
Se  as  3  500  m  one-sigma  and  2.300  m  SEP  for  time  errors  on  the  order  of 
1  seTone  sfgma  Discussions  with  systems  experts  suggest  that  time  errors  wfl 
be  on  ^  order  of  0.3  millisec.  Under  these  conditions 

on  the  order  of  1  m  one-sigma  for  Case  1  and  0.8  m  SEP  for  Case  2.  Clearly 
time  error  on  the  order  of  less  than  a  millisec  will  have  no  appreciable  effect  o 

track  accuracy. 
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Velocity  magnitude = 200  kt  Acceleration  magnitude = 3.5  g 


.^-"usec 


Velocity  magnitude =1,000  kt  Acceleration  magnitude = 9  g 


Velocity  magnitude = 1 ,600  kt  Acceleration  magnitude = 9  g 


4.  CONCLUSIONS 


4.1  Findings 

This  report  provides  the  details  and  the  results  of  an  analysis  technique  to  un- 
quantify  how  errors  in  navigation,  sensor,  and  time  affect  the 
SIAP  The  analysis  approach  derives  the  equation  for  the  location  of  an  aer 
s^ce  oblert  in  the  WGS  84  coordinate  frame  and  determines  how  errors  in 
measurements  and  alignments  affect  the  estonate  of  track  posmon.  This  r 
necessarv  since  meaningful  information  can  be  exchanged  between  the  snip 
and  othCT  units  in  a  distributed  system  only  if  the  location  and  orientation  o 
the  various  coordinate  frames  are  known  relative  to  each  other. 

An  examination  of  the  track  equation  identifies  three  groups  of  error  contribu- 

tors: 

•  Navigation  position  and  misalignments 

•  Sensor  measurement  and  misalignments 

•  Time 

After  the  error  mechanisms  are  identified,  a  first-order  linear  perturbation 
"«ed  to  develop  the  relationship  between  track  error  and  the 

errors  identified  above.  The  result  of  the  analysis  is  the  sens 

curaev  to  the  navigation,  sensor,  and  time  errors.  Once  the  sensitivities  are 

calculated,  the  error  budget  and  parametric  sensitivities  are  determine  . 

A  nominal  error  budget  for  the  navigation 

Frnm  this  the  track  accuracy  error  budget  is  developed.  This  budget  snows  a 
m  OToVconWbutlon  between  the  navigation  and  sensor  subsystems. 

“e  time  error  is  a  relatively  small  contributor,  so  that  m^“e  ^ 

from  the  nominal  value  will  result  in  small  decreases  m  track  error.  An  exam 

nation  of  the  navigation  subsystem  shows 

contributions  from  the  navigation  misalignments  and  the  lever  aim 
very  small.  The  sensor  subsystem  error  budget  is  almost  equally  dmd 
tween  the  sensing  errors  and  the  sensor  misalignment  errors. 

The  sensitivity  of  track  error  to  time  error  depends  on  the  aerospace  object  ve- 

todty  and  acceleration.  The  sensitivity  is  found  to  eq^  toe  rector  of  th 

velocitv  and  the  product  of  the  acceleration  and  the  known  time  latency  ot  tn 
system  For  a  high-speed  combat  aircraft  and  a  network  latency  of  5  sec 
SeTrror  is  maSified  by  a  factor  of  956  m/sec.  Thus,  a  time  error  of  !  sec  is 
i  t  a  ntn  a  track  error  of  956  m.  The  known  time  latency  of  the  system 
“"of  increasing  or  decreasing  toe  contribution  of  the  object's  aecel- 
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eration  to  the  sensitivity.  For  example,  reducing  the  time  latency  by  a  factor  of 
two  can  reduce  the  track  error  to  735  m. 

The  parametric  analysis  shows  that  worst  case  conditions  of  1  sec  time  error 

md  30  sec  network  latency  lead  to  track  errors  of  aPPro™"®‘e^°^™  £  er. 
high-speed  aircraft.  Time  errors  on  the  order  of  only  a  milhsec  lead 
rors  of  less  than  4  m  for  these  worst  case  conditions. 


4.2  Future  Work 

Although  the  analysts  assumes  a  shipbome  sensor  system 
eral  and  can  be  applied  to  all  systems  that  contribute  to  the  SIAP.  We  an ^ cur 
rently  developing  error  budgets  for  land  and  airborne  sys  ems.  n  ' 

are  emlnding  the  list  of  error  contributors.  For  example,  the  analysis  of  the 
time  en-or  In  lection  2.8  develops  error  expressions  for  extrapolation  and  tr 

We  are  investigating  how  these  errors  can  be  related  to 

tracking  errors. 

One  obiective  of  the  error  budget  analysis  is  to  determine  how  the  individual 
errors  ^Tect  the  SIAP  attributes.  We  are  investigating  how  to  “dorp°”‘  . 
error  analysis  approach  into  the  correlation/decorrelation  algorithms  to  de 
mine  sensitivities  of  the  attributes  to  the  individual  error  contributors. 

We  also  plan  to  incorporate  the  Common  Reference  Scenarios  (CRS1 1  into  this 
lalvS^ttiat  comparisons  can  be  made  to  other  simulation  results  based  on 
the  CRS.  Using  the  CRS  will  also  provide  insight  into  how  track  accuracy 
changes  as  the  threat  evolves  in  time. 

We  are  also  using  the  error  model  and  the  results  of  the  analysis  runs  to  de¬ 
velop  test  cases  to  assist  with  the  IABM  verification  and  validation  effort.  Th 
error  models  provide  the  mechanism  to  introduce  the  correct  level  of  bias 
each  system,  mid  the  results  of  the  error  model  analysis  provide  examples  of 
track  accuracy  that  the  IABM  can  be  checked  against. 
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